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and Scheffe's multiple contrasts to demonstrate how the TI 59 may be used 
to facilitate statistical inferences. The fitness data are also used 
for this purpose. The intent throughout the thesis is to exemplify the 
capabilities of the TI 59 as a viable, real world analytical tool rather 
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feet eODUeCTION 


The intent of this thesis is do demonstrate through three compre- 
hensive examples, the tremendous capabilities of the TI 59 orogrammable 
hand-held calculator. One of the examples is a probability application 
while the other two entail use of the TI 59 in statistical inference 
and data analysis. 

The example chosen to illustrate an application to probability 
theory 1s a combat simulation model. The model involves use of the 
Monte Carlo technique to simulate stochastically detection, identifi- 
cation, and engagement of a cruise missile by an Improved Hawk Air 
Defense Battery. Chapter III discusses the combat model and the TI 59 
Simulation in detail. Chapter II addresses briefly the Monte Carlo 
technique, combat modeling in general, and the TI 59 random number 
generator in order that the reader may better understand the combat 
model discussed in Chapter III. 

Two examples are provided for statistical applications of the 
calculator. These involve the analysis of data gathered by the authors 
in a physical fitness experiment. The experiment, discussed in detail 
in Chapter IV, involved the testing of thirty-six male subjects who were 
divided into six categories based on their training routine, to deter- 
mine whether the subjects’ training program did in fact influence their 
physical fitness. Chapter IV describes the fitness experiment in terms 
of its scope, design and scoring methodology. Chapter IV also discusses 
a TI 59 program which computes measures of central tendency and spread 


and then illustrates the use of the program with the fitness data. 
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Chapter V describes a TI 59 program for analysis of variance and then 
demonstrates how the program may be used with fitness data to 
make statistical inferences. 

Throughout the thesis, it is assumed that the reader is generally 
familiar with programming techniques for the TI 59 hand-held calculator. 
Subroutines, labels, flags, data registers, and program steps are 
discussed in each of the ensuing chapters where the intent is to 
illustrate how the features of the calculator may be exploited to 
Facilitate statistical analysis or simulation. Reference 12 discusses 
programming techniques for the TI 59. The TI 59 has one particular 
feature which makes it much more than a calculator. Specifically the 
capacity to use subroutines provides a analytical tool more like a 
minicomputer than a calculator. The three programs discussed in the 
succeeding chapters use subroutines extensively to illustrate this 
powerful capability. 

Finally, while a few of the referenced tables and charts of this 
thesis are positioned close to comments discussing their purpose, 
most are to be found at the conclusion of each chapter or in the 


appendices. 





Mit PROBABIUTDY 


The intention of the authors was to begin this chapter discussion 
with a definition of probability theory, that branch of mathematics 
generally believed to have been founded by a Swiss mathematician named 
Jacques Bernoulli. However, research has revealed that there is some 
discussion as to the true meaning of probability theory and that among 
mathematicians there appear to be those who view probability as a state 
of the universe while others consider probability a state of belief. 

To compound this situation furthermore, there appear to be differing 
definitions of probability within each group. Indeed, all attempts to 
define probability directly have failed to meet with success. Instead, 
probability has been axiomitized, much like geometry, so that a set of 
consistency rules or axioms established by A. Kolmogorov are now 
generally accepted by the scientific community. These axioms allow a 
great deal of freedom in the assignment of probabilities for any par- 
ticular model and at the same time force any such assignment to be 
consistent with any other. Moreover, the theorems of that theory then 
become universally true statements for any such assignment. In this 
system, events are defined as sets in a specified sample space. With 
those guidelines as a background, probability theory can be used to make 
intelligent predictions and decisions if we know what events are possible 
and how probable are the various events. After a little thought it be- 
comes immediately apparent that the immense power of such a tool as 


probability theory is limited in use only by one's imagination and 


10 





ingenuity. This research is an effort to use probability theory in the 
construction of a probabilistic combat simulation on the Texas Instruments 
programmable 59 calculator (TI 59). Because the simulation developed 
includes a number of the many chance elements involved in most combat 
Situations, a discussion of the Monte Carlo technique and random number 
generation on the TI 59 follows. A brief disucssion of combat model 
Simulations concludes this chapter. 


A. MONTE CARLO TECHNIQUE 


sf 
a 


Systems that exhibit stochastic elements in their behavior can be 
Simulated with the aid of the technique called Monte Carlo (named after 
the famous gambling resort town of Monaco). This technique involves 
Sampling from those known probability distributions that represent each 
of the actual chance processes included in the system under study Es 3]. 
By completing a system simulation run many times while keeping the non- 
Stochastic inputs constant but allowing the chance elements to fluctuate 
according to their known probability distributions, a statistical 
average for run results can be determined. 
Turban and Meredith free. 3| have listed the steps necessary in 
building a Monte Carlo simulation as follows: 
"lL. Describe the system and obtain the probability distributions 
of the relevant probabilistic elements of the system. 
2. Define the appropriate measure(s) of performance. 
Beeconstructe Cumulative probability distributions for each of 
the stochastic elements. 
4. Assign representative numbers in correspondence with the 


cumulative probability distributions. 
5. Generate a random number for each of the independent stochastic 


elements and .. . (determine) the measure of system performance. 
6. Repeat step five until the measure of system performance 
Stabilizes." 


ugk 





Thus the distinguishing feature of the Monte Carlo method is the re- 
petitive execution of an established experiment or simulation involving 
randomness. 

While electronic digital computers themselves are not necessary for 
the execution of simulations, they do offer tremendous speed and con- 
Sistency of conditions for such models. Thus the computer is ideally 
suited to perform the large number of repetitions required by Monte 
Carlo but the matter of randomness presents a problem. For the Monte 
Carlo technique described above the necessity of a truly random number 
is essential. However Kovach [net 6| notes that: 


"Strictly speaking, the random number exists only as the 
result of a random process." 


While computers, to include the TI 59, do possess the capacity to con- 
tinuously generate random numbers as they are needed, these numbers are 
subject to the limitations of the computer and are not truly the result 


of a random process and hence are often described as pseudo-random. 


ee TI 59 RANDOM NUMBER GENERATOR 
R.F. Barton Res. ‘| describes simulation as follows: 


"Simulation is simply the dynamic execution or manipdu- 
lation of a model of an object system for some purpose. 


Simulation is a case-by-case method for studying 
object systems. Each case might be either a single trial 
or an entire run. In either view, outputs may differ 
meal tO trial and run to run.” 
The object system is that system under study in the simulation. 
The TI 59 in its capacity as an electronic computer provides the 


user with the means of developing and executing stochastic and nonstochastic 


Simulations. 
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Barton continues Ref. 4| : 

" A stochastic simulation is one in which differing outputs 
trial to trial can be obtained without changing the inputs 
(ignoring random numbers as inputs). Specifically, this 
means that identical parameters, starting conditions, and in- 
put time path values produce varying outputs trial to trial 
aie on CO run. 

A nonstochastic simulation is one in which the inputs or 
the model must be changed to obtain changed outputs. This 
means that identical model operations, parameters, starting 
conditions, and input time path values will produce identical 
SUeputs run to run. " 

There are inputs common to both of these simulation types. However, 
as alluded to above, there are also special inputs that are needed to 
represent the chance processes or stochastic events found only within 
a stochastic simulation. These special inputs are random numbers. 

The characteristic of random numbers that makes them different from 
all other numbers is the fact that the knowledge of any future random 
number cannot be enhanced by the knowledge of any past, present, or 
other future random number. 

The TI 59 with its master library module solid state software pro- 
gram ML-15, a random number generator, can generate sequences of uniformly 
or normally distributed random numbers independent of a simulation pro- 
gram or within such a program. 

Kovach states [Ree. é| : 

"(Random) numbers generated by the computer are sometimes- 

called pseudo-random because they are subject to the 

limitations of the computer. In a list of truly random 
numbers, for example, one would expect to find numbers 
containing more digits than can be obtained in a computer.” 


Random numbers produced by the TI 59 ML-15 program are generated by 


a mathematical formula. Given an initial seed number by the user, this 
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program will always produce the same list of pseudo-random numbers. 
Thus if repeatedly initialized with the same seed number the forth- 
coming random numbers would be known and randomness would not exist. 
That is, every future random number could be predicted. Hence, the 
randomness of the numbers produced by the ML-15 program are as dependent 
upon the user as the mathematical formula of the program itself. It is 
therefore incumbent upon the routine user of the ML-15 program to vary 
the seed number used within denoted limits to insure genuine pseudo- 
random numbers. 

The TI 59 ML-15 random number generator program is listed in 
Table 2-1. User instructions for the ML-15 program Ref. 1] are 
listed in Table 2-2. Data register contents are listed in Table 2-3. 

While the program does provide the option of generating uniform 
numbers for ranges other than 0-1 and also provides statistical data 
for the random numbers generated, only the generation of the uniform, 
range O-l, and the normal random numbers will be discussed further. 

1. Uniform Random Number Generator 

A flowchart of the uniform random number generator, range O-l, 
1s displayed in figure 2-1. Program steps 000 through 054 contain the 
following mathematical formula, called the linear congruential 
method Ree. 1 , for the generation of these numbers. (Throughout 
this thesis an asterisk is used to indicate multiplication.) 
( (24298 _ SEED - 99991) — 199017 STO 07) 

A work value is the result of the above operation. This result remains 


in the display register. The value 199017 is stored in data register 


O7, an ML-15 work register. 
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NUMBER GENERATED 
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New SEED 
GENERATED AND 
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STOP 


TI 59 Uniform Random Number Generator 


PIGURE 2-1 


ile, 
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Galculation continues. 

(INV INT i RCL) O77) STO 69 
The integer portion of the number resulting from the previous operation 
is discarded, then the remaining fractional portion is multiplied by 
199017 which was stored in data register 07. This product then is 
stored in data register 09 and becomes the seed for the next random 
number calculation. Calculation continues. 

((RCL O09 — RCL 07) * 5S INV 2nd log) 
Now the new seed is divided by the number 199017 which was stored in 
data register 07 during the first operation. This quotient is then 
multiplied by the common antilogarithm of 5 to complete the step. 
Sareulation continues. 

(INT — 5 INV log) 
The fractional portion of the previous numerical operation is discarded 
and the result divided by the common antilogarithm of 5. This quotient 
is then displayed as the uniform random number, range O-l. 

2. Normal Random Number Generator 

A flowchart of the normal random number generator is displayed 
in figure 2-2. Program steps 069 through 135 contain the following 
mathematical formula, called the direct method Res. Ly for the 
generation of these numbers. 

(seed STO 09) 

(mean STO 10) 

(standard deviation STO 11) 
Following program initialization, three data values are entered and 


Stored in data registers 09, 10, 11. The seed number, the desired 


normal distribution mean (mu) and the desired normal standard deviation 
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(sigma) are stored in these registers as listed respectively. The seed 
value is limited as noted in the TI user instructions, table 2-2. 

(SBR DMS STO O08 SER DMS) 
Initially the program calls the previously discussed uniform random 
number generator to produce such a number and then stores it in data 
register 08. The uniform generator is called again to produce another 
uniform random number which remains in the display register for manipu- 
lation and is denoted RN. Calculation continues. 

Giaie * 32 x “i } “EOS) = ((RCL O8 LnX) * (-2) 
The uniform random number in the display is next multiplied by two pi. 
The cosine of this product is then calculated. The resulting value is 
multiplied by the product of the natural logarithm of the first uniform 
random number (generated early and stored in data register 08) multiplied 
by negative two. This product remains in the display register for 
manipulation. Calculation continues. 

(GV * RCL 11) + °9RCL 10 
The squareroot of the previous operation end value is multiplied by the 
desired standard deviation. Finally, this product is added to the 
desired mean, resulting in the generated normal random number. Seed 
manipulation for the generation of successive normal random numbers is 
completed during the SBR D.MS portion of the normal generation program. 

Ber obR cna. OMS 

The TI 59 ML 15 program does compile statistical data to allow 

computation of the mean and standard deviation of the pseudo-random 


numbers generated when using the normal distribution routine and the 
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uniform distribution routine over ranges other than 0-1. However, these 
data are not compiled when SBR 2nd D.MS sequence is executed to pro- 
duce uniformly distributed numbers over the O-l range. Data registers 
one through eleven are used by the ML-15 program to compile and compute 
these statistics. Hence, if this program is called to produce normal 
random numbers within a larger program, such as a simulation, the use 

of these eleven data registers must be forgone. Yet, if the ML-15 voro- 
gram is called only to produce uniform random numbers over the range 
O-1, only data registers seven and nine are used by ML-15, freeing 

nine registers for other use. This aspect of the TI 59 ML-15 program 
must be carefully considered when utilizing it is as a subroutine within 


another program. 


C. COMBAT MODELS 

Today there are considered to be three types of combat models in 
use; War games, pure simulations and analytical models. War games are 
models and we ate [Ref 2| : 


. « « in which individuals simulating decision makers 

in real life use their judgement to verform the 

decision functions in the model." 

A war game may include automation to assist in the processing of data 
and the generation of random numbers to determine the outcome of certain 
chance events. A war game may also be a player-assisted simulation 
where players provide input to a computer model based on output (read- 
outs) during a simulated battle. In comparison with the other models, 
war games appear to be more realistic, involve greater player interac- 


tion, are less automated, require much more time to run, more resources 


and involve a smaller degree of abstraction. 
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Simulation combat models are models [Ree 2| : 


". . . which run completely without human intervention. 


In this type of model events in the different combat 

processes are based on predetermined rules which are 

programmed into the automated evaluation procedure." 

Combat models of this type generally contain a significant number of the 
important stochastic elements of combat in an attempt to simulate real 
battle. These models use probability distributions for the many chance 
input variables and produce probability distributions as results. 

They utilize the Monte Carlo technique, repeatedly sampling all input 
distributions in the programmed sequence to produce a distribution of 
probable battle results for each set of input data. 

Analytical models are models Es 2| ; 

" . . . Comprised of sets of mathematical equations as 

models of all the basic events and activities in the 

process being described and an overall assumed 

mathematical structure of the process into which. the 

event or activity descriptions are intecrated. " 
While analytical models are the most time efficient they are also the 
most abstract and difficult to understand. As with the pure simulation, 
there is no human intervention when an analytical model is used. 

All three models represent abstractions of the real world. The 
models themselves can be observed more conveniently than the real world 
and theories’ about the real world can be developed by studying the re- 
sults of these models. Subsequently, these theories can be used to 
make predictions about real world events. 


Each model type has strengths and limitations, some noted above and 


others listed in table 2-4 Ee 2| A 
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1. Pure Combat Simulations 

Pure combat simulations are normally viewed as production tools, 
using Monte Carlo techniques to obtain results enabling the prediction 
of future system performance. But because the real world is so complex 
and interactive, attempts to model every detail of a large system in a 
pure simulation and to include every element that may influence the 
system can result in simulations so large and so complex that they 
are understood only by their developers and not by other users or 
decision makers. To avoid this complex dilemma, analytical models 
can be used to represent elements of the system being modeled instead 
of simulating the element itself and its inherent stochastic processes 
with every trial. This technique has been followed to some degree in 
the pure simulation model presented in Chapter III where the calculation 
of detection probability is an analytical model with results based 
Sameange: range. 

It should be noted that few, if any, simulation models ever com- 
pletely include all those elements and events that affect the system(s) 
under study. Reference 2 points out that: 

", . . a model is always incomplete, with only those 

aspects represented that we believe we know well 

enough to model and that we consider important in 

the issues to be examined with the model. Obviously, 

models tend to be as simple and concise as our 

knowledge of the activity warrants." 
This is reflected in the model presented in Chapter III. While all the 
factors affecting system performance have not been directly simulated 


they are included either as analytical models or as given in the 


scenario. 
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Finally, a point to be stressed is that simulations need not be 
Yarge to be useful, nor require the use of a large electronic digital 
computer to be credible. Using large computers for large problems and 
small computers for small problems is a rule of thumb that may over- 
state the case but certainly does not exaggerate it. Use of the TI 59 
aS a computer to tackle the problem set forth in Chapter III is an 
example of matching the problem to the computing power required. It 


1s also an excellent example of the computing power of the TI 59. 
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TI 59 PROGRAM ML-15 
USER INSTRUCTIONS 


(MASTER LIBRARY MODULE) 


PROCEDURE ENTER DISPLAY 


SELECT PROGRAM 2d PGM 15 














INITIALIZE Zee Q 





ENTER RANDOM NUMBER SEED SEED E SEED 





(OS SEED = 199017) 


FOR NORMAL DISTRIBUTION 
ENTER DESIRED MEAN Xx A 
ENTER DESIRED STANDARD 

DEVIATION on B 
GENERATE RANDOM NUMBER 

(REPEAT AS NEEDED) 2a 1C 
FOR UNIFORM DISTRIBUTION 

RANGE (0, 1) 
GENERATE RANDOM NUMBER 

(REPEAT AS NEEDED) SBR 2d D.MS 


TABLE 2>2 
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TI 59 PROGRAM ML-15 


DATA REGISTER CONTENTS 


BEGUESTERS: 
Roo = 
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Roo ya 


o2 «ON 
Rog xX? 
Rog LX 


Rog 2 XY 


Be WORK REGISTER 
an WORK REGISTER 
Rog SEED 
Rig A,% 
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War Games 


Intelligent play of decision maker. 


Intelligent use of intelligence. 


Can plan engagement and moves in 
advance. 


Adapts maneuver to situation. 

Can play many tactical situations 
(employment, penetration, etc.). 
Insight is gained by understanding 
the particular rationale used in 
the decision process in single 
Situations analyZed. 

Controller determines existence 

of engagement and pace of play. 


Very slow and costly in resources. 


Very few situations can be 
examined. 
Greatest visibility for the user. 


Can include direct involvement by 
user. 


Simulations 


Stvirzed decision routines, 
usually fixed throughout game. 


Very limited use of intelligence. 


Very limited planning horizon. 


Very limited adapted maneuver 
routines. 


Usually stylized maneuver, limited 
change in formations. 


Insight 1s gained by repeating 
the analysis in many situations 
uSing different values for key 
parameters. 


Predetermined scenarios and 
engagement rules -- combat very 
intense. 


Faster to run after completely 
developed. 


Many situations can be played 
and the sensitivity of key 
variables can be tested. 


Reference 1, FM 44-90, Headquarters Department of the Army, 1977 


TABLE 2-4 
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[lie COMBAT SIMULATION USING TI-59 


This is a probabilistic duel simulation model, a pure simulation of 
a combat air battle, designed to reflect the characteristics of the 
Improved HAWK air defense artillery system (battery) in the manual 
mode under attack by a single cruise missile of sustained altitude, 
speed and direction. 

The program scenario and engagement rules are predetermined with 
no user input once the simulation run has begun. Insight may be gained 
and the sensitivity of key variables tested dependent upon the use of 
different values entered by the user for these variables during program 
initialization. (Variables listed-under E below.) 

The model provides IHAWK system status, target engagement events and 
battle results as they are determined/occur. Only two battle results 
are possible: a "KILL" of the cruise missile or a unit “PENETRATION" 
by the cruise missile. 

Given the operational ready rates of the major subsystems of the 


IHAWK system, the P (probability of single shot kill) and the mean 


ssk 
and standard deviation of lock-on-to-target times, the model samples 

from the uniform and normal distribution to determine system status, 

IHAWK missile kill or no kill and lock-on times. Target detection is 
modeled as a function of target range and is represented as a linear 

relationship in the simulation. 


This simulation was developed to exhibit the computing power of the 


TI 59 and to determine if one parameter under the control of the IHAWK 
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battery commander could significantly affect air battle results. This 
parameter was the tracking radar "lock-on-to-target"™ time which is a 
function of operator training given (1) a manual mode operation, and 
(2) perfect equipment. The sensitivity of battle results to varying 


lock-on times is listed under F below. 


A. IHAWK SYSTEM 
To be effective, an air defense system must be able to detect, 

identify, engage and destroy an airborne target. The IHAWK system can 
engage and destroy a full spectrum of threat aircraft and missiles operat- 
ing throughout a wide range of tactical speeds and altitudes. It can 
engage a multiple target threat as well as single targets. The system 
is effective from ground level up to altitudes of about 48000 feet and 
out to ranges of about 40 kilometers. The system can operate at night, 
under all conditions of weather and reduced visibility. It can function 
effectively in an ECM (electronic countermeasures) environment and is 
mobile using organic unit vehicles or helicopters [Ree 1| . However, 
some adverse weather and heavy ECM may diminish some system capabilities. 
This simulation does not model: (1) weather 

(2) detection ECM 

(3) viability 

(4) system mobility 

1. Detection 
Target detection is accomplished by either the improved pulse 

acquisition radar (IPAR) or the improved continuous wave acquisition 
radar (ICWAR), or both. The IPAR can detect low to medium altitude 


targets out to ranges in excess of 100 kilometers while the ICWAR can 


detect targets at very low altitudes with ranges in excess of 60 








a> Detect 
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kilometers. Operating on the doppler principle, the ICWAR sees only 
very low moving objects and not stationary objects on the ground. 
Detection may be reduced if terrain features such as hills, trees and 
buildings mask the radar's view of the approaching target. Even with 
completely level terrain, the earth's curvature causes a reduction in 
the detection and tracking ranges of the system against very low-altitude 
targets. Also, evasive maneuvers by threat targets may reduce the de- 
tection and tracking ranges and increase system response time, thereby 
reducing the effective or intercept range. This simulation does not 
model: (1) terrain features 

(2) curvature of the earth 

(3) evasive action by targets 

(4) pulse detection or continuous wave detection per se 
This model assumes: 

(1) clear weather 

(2) no detection electronic countermeasures 

(3) line of sight (LOS) exists between radar and target 

(4) flat desert terrain 

(5) nonmaneuvering target 

(6) only one attacking target exists 

(7) target is a cruise missile of constant speed 

and constant altitude 

(8) detection is a function of target range 
The probability of detection is modeled during each sweep of the radar 
as a linear function of target range from the battery as follows: 
Hemeene IPAR the Popp = |(-.25 + 65). * Target Range| => S20 
For the ICWAR the Ppp = [i-.5 = 65) * Target Range| + 1.0 
The probability of detection is calculated every three seconds of 
Simulated time. This is based on the radar rotation rate of 20 


revolutions per minute. That is, every three seconds each radar 


takes a 360 degree glimpse of the horizon. The radars are slaved to 
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each other and rotate in synchronization. Additionally, the IPAR is 
modeled to detect only targets from 5000 to 40000 feet in altitude while 
the ICWAR detects targets from 1 to 8000 feet in altitude. Thus the 
battery's very low and low to medium detection capability is dependent 
upon the operational status of these radars as noted. Targets above 
40,000 feet cannot be detected in this simulation. The operational 
ready rates of these two radars has been arbitrarily set at .65 (ICWAR) 
and .95 (IPAR). 
eee Laentification 

Identification of any potential target is accomplished by means 
of the identification, friend or foe (IFF), equipment of the IHAWK 
system and/or other established hostile criteria. If the target can- 
not be positively identified in this simulation because of a non- 
operational IFF, the target speed and altitude is checked to determine 
target status (foe or not foe). That is, if the IFF is nonoperational 
and if the target is below 5000 feet altitude and greater than 550 KMPH 
in speed, it is identified as a foe; otherwise, it is not a foe. 


This model assumes: (1) once identified as a friend, always a friend, 
(2) once identified as a foe, always a foe. 


The operational ready rate of the IFF has been arbitrarily set at 

95 percent. This model does play IFF accuracy to the degree that an 
operational IFF will be in error two percent of the time. That is, a 
foe will be shown to be a friend two percent of the time. This model 


assumes operational IFF accuracy to be 98 percent. 
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3. Engagement/Destruction 

The IHAWK battery has two firing sections, each of which 
contains a target tracking radar which is called an improved high powered 
illuminating radar (IHIPIR), and three launchers with three missiles 
each. If a potential target is determined to be a foe, it is assigned 
to one or both of these firing sections. The tracking radar of these 
sections, under the control of a fire control operator (enlisted 
personnel), attempts to lock-on to the approaching hostile target. The 
operator directs the automatic box search of the radar in the azimuth and 
expected elevation of the target in attempting this target lock. The 
operator's ability to achieve a target lock is a function of his train- 
ing given the condition mentioned above. The time that elapses in 
attempting this target lock is extremely important. Target engagement 
cannot continue without target lock and the longer it takes to achieve 
target lock the closer the target moves toward the battery, reducing the 
intercept range. Tracking radar lock-on times are assumed to be normally 
distributed. After target lock has been achieved and the target is in 
range, one or two missiles are fired on order from the unit tactical 
control officer in a battery control van. Engagement is continued 
until the target is destroyed or until engagement is no longer possible. 
This simulation models each firing section as an entity. After a 
target has been declared a foe and assigned to one or both firing 
sections, this model simulates the target lock-on time by utilizing the 
TI 59 normal random number generator based on a normal mu and sigma 


input by the user during initialization. 
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The range at intercept is determined by the range at which detection and tracking 
(lock-on by the HIPIR) occur and on system response time. 
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In this model, targets are engaged that are: 
(1) declared to be foe 
(2) less than 40 KM from the battery 
(3) greater than 8 KM from the battery 
\ 
This model assumes: 
(1) two independent firing sections 
(2) salvo fire occurs if both sections are operational and 
shoot-look-shoot if only one section is operational 
(3) firing continues until kill or penetration 
(4) penetration means that the target is 8 KM or less from 
the battery 
(5) lock-on-to-target time is a function of operator 


training and is normally distributed 
(6) P = .75 (arbitrarily set) 


ssk 
The operational ready rate of each firing section has been arbitrarily 
set at 75 percent. 
ee berget 

The target for this model is assumed to be a hostile cruise 
missile that flies straight in toward the battery at a constant speed 
and altitude as established by the user during initialization. The 
initial range of the target is also a user input. The lethality of 
the missile warhead is assumed to be such that any successful penetra- 
tion by the missile to within 8 KM or less of the battery before 
destruction is considered a total penetration of the battery defended 
area. Therefore, the target must be destroyed before 8 KM to score a 
kill. Additionally, a target will not be engaged after detection until 
it is less than 40 KM from the battery and no further missiles will be 
fired at the target once it is within 8 KM of the battery. The target 
speed has a lower bound of 100 KMPH but no upper bound. Only targets 
between 1 and 40000 feet in altitude can be detected and are thus the 


altitude bounds. Finally, all targets are hostile and will be engaged 
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mmbess erroneously identified as friendly or not foe, resulting in a 
free penetration. 
>. Time 
This simulation is a time step model, updating all battery 
events and functions every three seconds of simulated time. This three- 
second interval stems from the rotation rate of the detection radars, 
20 revolutions per minute or one complete rotation (scan of the horizon) 


every three seconds. 


B. MACRO FLOW CHART 

The enclosed macro flow chart, figure 3-4, depicts the general flow 
of the simulation logic from start to either penetration or kill. 

First the model determines if a detection capability exists. This 
could be one or both of the detection radars. Using the internal random 
number generator of the TI 59 for a O-l uniform distribution, two 
random numbers are drawn and compared with the detection radar opera- 
tional ready rates. If the random number is less than the rate, the 
radar is operational; otherwise, it is nonoperational. If no detection 
capability exists the simulation is terminated by a penetration of the 
defended area by the target. 

Detection of the target is based partly on the formulas set forth 
above and results of the O-1l range uniform random number generator. 

The probability of detection is based on the range of the approaching 
target and is recalculated every three seconds. The probability of 
detection for each radar is compared with a generated random number 


between 0 and l. If the random number is less than the probability 
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the target is detected; otherwise, it is not. As the target moves 
closer to the unit, the probability of detection increases linearly. 
However, detection is not based on range alone. The target must be 
within the detection altitude capability of the operational system. 

For instance, if the battery's detection capability rests solely on the 
IPAR because of a nonoperational ICWAR and if the approaching target is 
at an altitude of less than 5000 feet, detection will never occur and 

a penetration will result. 

Identification of the target is determined by either the IFF or 
a combination speed and altitude envelope if the IFF is nonoperational. 
Operational status of the IFF is determined by comparing, again, a 
generated random number from a O-1 range uniform distribution against 
the operational ready rate of the IFF. If the random number is less 
than the ready rate, the IFF is operational; otherwise, it is non- 
operational. If the IFF is nonoperational and if the unit has a 
tracking capability, the target speed and altitude can be checked 
against an established hostile criteria. If the target is found to be 
below 5000 feet in elevation and above 550 KMPH in speed, it is 
designated a foe; otherwise, it is considered to not be a foe. Only 
targets positively identified as foes are engaged. 

IFF positive identification of the hostile cruise missile as a foe 
is determined by comparing a random number generated from a O-l range 
hierorm distribution against the IFF accuracy rate of 98 percent. If 
the random number is less than the accuracy rate, the cruise missile 
is correctly identified as a foe; otherwise, it 1s erroneously classi- 


fied as a friend resulting in a penetration of the defended area. 
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The battery's engagement capability lies with its two firing 
sections, each tracking and firing on approaching hostile targets. 
Again, two random numbers drawn from the 0-1 uniform distribution are 
compared to the operational ready rates of the firing sections to de- 
termine the status as in the previous subsystems examples. If an 
engagement capability is determined to exist, the process continues on 
to direct missile firings at the target following target lock by the 
tracking radar(s). As mentioned earlier, this lock-on time is a 
function of operator training in the manual mode and is normally 
distributed. For highly trained operators the mean is assumed to be 
ten seconds with a standard deviation of five seconds. The lock-on 
time for each section is determined by two random numbers generated by 
the TI 59 random number generator from a normal distribution with mean 
and standard deviation determined by the user during program initializa- 
tion. (Any random number less than zero is discarded and another 
generated to avoid negative times.) The length of the lock-on period 
directly affects the resultant target intercept range. If the lock-on 
time is sufficiently long the missile firing is delayed and the proba- 
bility of a penetration is likely. (The determination of air battle 
results to varying lock-on times provided the basis for the development 
of this simulation, though other variables of the model can easily be 
tested for outcome sensitivity.) Again, targets are not engaged until 
less than 40 KM from the battery and no missiles are fired after the 
target is 8 KM or less from the battery. 


Missile effectiveness after firing is determined when the IHAWK 


missile range equals or exceeds the cruise missile range from the 
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battery. A random number from a O-1 range uniform distribution is 
generated for each missile fired and checked against the P,.,. 
femaom numbers less than the P..;, result in kilis while all others 


result in no kills. Engagement of the target continues until a kill 


Or penetration is registered. 


C. MICRO FLOWCHART 
The micro flowchart in appendix A depicts the detailed flow of 


processing throughout the simulation from start to kill or penetration. 


D. SUBROUTINES, LABELS, FLAGS, DATA REGISTERS AND PROGRAM MEMORY STEPS 
1. Subroutines and Labels 
This simulation uses 49 of the 72 labels available for pro- 
gramming on the TI 59. Of these 49, 14 are subroutines. The remaining 
labels are used to identify sections of the program and to direct 
action to these sections during simulation runs. A complete listing 
of all labels is displayed in table 3-1 with the subroutines marked 
by an asterisk. Comments on selected labels and a printout of the 
entire program is enclosed in appendix A. 
2. Flags 
Nine of ten available flags are used. As the IHAWK equipment 
and target friend/foe status is determined, this model uses TI 59 flags 
to maintain a record of the system and target status. These flags 
subsequently direct the flow of processing and determine actions to 
be taken within the simulation engagement. 
Flag 1 set means the target has been detected. Flag 2 set 


means the target is a friend. Flag 3 set means the ICWAR and the IPAR 





are operational. Flag 4 set means only the ICWAR is operational. 
Flag 5 set means only the IPAR is operational. Flag 6 set means the 
IFF is operational. Flag 7 set means that alfa firing section is 
operational and bravo firing section 1S nonoperational. Flag 8 set 
means that bravo firing section 1S operational and alfa firing section 
is nonoperational. Flag O set means that both firing sections are 
operational. 

See Data Registers 

The TI 59 memory storage area is initially partitioned to provide 
60 data storage registers and 480 program storage locations. However, 
the user can repartition the memory storage area to suit his particular 
programming needs. The IHAWK simulation requires exactly 800 program 
memory locations and 20 data storage registers. Within the TI 59 there 
are a total of 120 registers to be used for data storage and program 
locations. While each register can store only one datum point, each can 
store eight program instructions or steps. Thus 8 * 60 = 480 program 
locations which are initially available as mentioned above. Reparti- 
tloning the core 120 registers is done in increments of ten. Hence, to 
get the 800 program steps for the IHAWK simulation 100 core registers 
are needed. This leaves exactly the 20 needed for data storage. 

To partition the storage area, the number of sets of 10 data 
registers needed is entered and 2nd OP 17 pressed. Thus for the IHAWK 
simulation, twenty data registers are available after the initial repar- 
titioning by pressing 2 2nd OP 17. The registers and their contents are 


listed on the following page. 


2 





4. 


Target range. Entered by user, 


Rg, through R 
program. 


og 2re used by the TI 59 random number generator 


SEED for random number generator. Entered by user. 

Mean lock-on-to-target time. Entered by user. 

Standard deviation of lock-on-to-target time. Entered by user. 
THAWK missile range from battery. Initially zero. 

THAWK misSile range from battery. Initially zero. 

Probability of detection work register. 

Target speed. Entered by User. 

Target altitude. Entered by User. 

Target range work register. Not entered by user. 

Simulation trials or runs to be completed. Entered by user. 


Simulated time in seconds for each trial. Initially zero for 
each trial. 


Program Memory Steps 


There are 800 program steps available. All 800 program memory 


steps are used in this program. 
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Bee USER INSTRUCTIONS 

The enclosed user instructions, table 3-2, provide the necessary 
steps to initiate a sequence of simulation runs. The enclosed printout 
results, table 3-3, indicate the 22 possible print statements that 
may occur during the simulation. A sample of data input and simulation 
run results are displayed in tables 3-4 and 3-5. 

Two steps of the user instructions warrant further comment. 

moe oceo 2s Clear Data Registers. 

Instead of clearing all data registers the user may wish to 
clear selected registers when repeating simulation runs as in the case 
of sensitivity analysis work. In this instance the user may just clear 
registers Roy Rio: Riz: Rize Rigs Rio and enter the desired values. 
Rog need not be reentered as the program automatically changes the 
seed after each random number is generated. If the user does clear 
all data registers with 2d CMs the user must then enter an entirely 
pteerent Seed in Rog within the bounds noted. 

2. Step 9 Check Data Register Content 

This step is a quick safeguard for the user to ensure that the 
Simulation run is based on the correct parameter values. This step 
provides a complete listing of the 20 data registers with contents for 


review prior to the final user step. 


F. SENSITIVITY ANALYSIS RESULTS FOR LOCK-ON TIMES 
Four hundred simulation runs were made with four lock-on mean 
values: 10, 20, 30 and 40 seconds. In each case the lock-on time 


Seamearcd deviation was five seconds. For each simulation run the 
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target was initially set at a range of 50 KM in Roos target speed was 
1500 KMPH and target altitude was 7000 feet. An initial random number 
seed was entered for run number one but no further user seeds were 
provided, thus leaving seed manipulation to the program. 


@ae results listed below indicate that air battle results are in- 


deed sensitive to target lock-on times. 


MmeCK-ON TIME PENETRATION KILL 
mu sigma 

10 5 11% 89% 
20 5 14% 86% 
30 5 26% 7 42 
40 5 44% 56% 


The results indicate a Significant increase (12 percent) in defended 
area penetrations for a mu of 30 under the present scenario. This trend 
continues at an apparent exponential rate. With a mu of 40 seconds, 
defended area penetrations increase another 18 percent. Based on these 
results it appears advisable to maintain such a state of operator 
training that the mean target lock-on times be twenty seconds or less 
with as little deviation among the operators as possible. Furthermore, 
it seems that for the extra training assumed to be required to reduce 
mean lock-on times from 20 to 10 seconds there appears to be only a 
small marginal reward in the reduction of defended area penetrations 


(3 percent). 


G. RECOMMENDATIONS 
While the intent throughout the thesis is to exemplify the capabili- 
ties of the TI 59 as a viable, real world analytical tool, the results 


of the TI 59 simulation lend insight into an area that requires 
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7 


further investigation, that being IHAWK target lock-on times. While 
only the lock-on times themselves were varied for this simulation, other 
important Scenario parameters should be varied to acquire an improved 
understanding of how air battle results can be affected by lock-on times. 

Future enhancements of these results would include a significant 
increase in simulation runs for a wide variety of Scenario parameter 
settings. While this TI 59 model allows certain parameter variations 
during program initialization, other parameters such as acquisition 
radar altitude detection capabilities can be varied with only minor 
adjustments to the program. 

Regardless of whether future simulation studies are conducted using 
this TI 59 model or a facsimile on another computer, the results 


mipeve warrant further research in this area. 
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USER INSTRUCTIONS 


SiLEP PROCEDURE ENTER PRESS DISPLAY 


REPARTITION 2d OP 17 


800 PROGRAM MEMORY 
LOCATIONS 

20 DATA MEMORY 
REGISTERS 


Enter magnetic program 
cards, sides 1 through 
4 


CLEAR DATA REGISTERS 

RESET ALL FLAGS AND 

CLEAR ALL SUBROUTINE 

RETURN REGISTERS 

ENTER IHAWK SYSTEM DATA 
TRACKING RADAR 


“EOCK-ON" TIME SLO 0 
Sto yul 


ENTER CRUISE MISSILE GT@./77 LRN 
DATA 


RANGE (KM) RANGE 

10 < KM < 100 

SPEED SPEED 
ALTITUDE ALTITUDE 
(FT) 


ENTER SEED FOR RANDOM SEED 
NUMBER GENERATOR 
(OS Sa4e1929017) 


ENTER DESIRED NUMBER 
OF SIMULATION RUNS 





Table 3-2.1 


Di 





USER INSTRUCTIONS 
SEP PROCEDURE ENTER PRESS DISPEAY 


9 CHECK DATA REGISTER 
CONTENT 
10 START SIMULATION RUN (SEE 
RESULTS 
POSSIBLE) * 
* PC-100 C PRINTER REQUI 





Table 3-2.2 
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NUMBER PRINTED 


i: 


LO 


ina 


2 


14 


5 


17 


3 


2S 


24 


25 


RESULT PRINTOUTS 
MEANING 


CRUISE MISSILE HAS BEEN DETECTED AT (RANGE GIVEN IN 
KMs FROM UNIT) 


LOW ALTITUDE DETECTION RADAR (ICWAR) AND MEDIUM 
ALTITUDE DETECTION RADAR (IPAR) ARE BOTH OPERATIONAL 


ONLY ICWAR IS OPERATIONAL 
ONLY IPAR IS OPERATIONAL 
IDENTIFICATION FRIEND-OR~-FOE (IFF) IS OPERATIONAL 


ALFA FIRING SECTION IS OPERATIONAL, BRAVO FIRING 
SECTION IS NONOPERATIONAL 


BRAVO FIRING SECTION IS OPERATIONAL, ALFA FIRING 
SECTION IS NONOPERATIONAL 


BOTH FIRING SECTIONS ARE OPERATIONAL 
CRUISE MISSILE IS IDENTIFIED AS A FOE 


BRAVO FIRING SECTION IS FIRING ONE MISSILE AT A TARGET 
<40 KM FROM THE BATTERY, BUT GREATER THAN 8 KM 


ALFA FIRING SECTION IS FIRING ONE MISSILE AT A TARGET 
<40 KM FROM THE BATTERY, BUT GREATER THAN 3 KM 


CRUISE MISSILE IS ERRONEOUSLY IDENTIFIED AS A FRIEND 
BY IEE 


CRUISE MISSILE IDENTIFIED AS NOT FOE BY SPEED AND 
ALTITUDE CRITERIA, IFF IS NONOPERATION 


ALFA FIRING SECTION IS OPERATIONAL 

ALFA FIRING SECTION IS NONOPERATIONAL 

BATTERY IS NONOPERATIONAL, NO DETECTION CAPABILITY 
BATTERY IS NONOPERATIONAL, NO FIRING CAPABILITY 


ALFA SECTION MISSILE "NO KILL" FOLLOWED BY RANGE (KM) 
OF APPROACHING TARGET 


Table 3-3.1 


a, 





NUMBER PRINTED 


26 


66 


RESULT PRINTOUTS 
MEANING 


BRAVO SECTION MISSILE "NO KILL" FOLLOWED BY RANGE (KM) 
OF APPROACHING TARGET 


IFF IS NONOPERATIONAL 


DARGE Lo KEonS JARE SPELLED OUT, I.E. “KILL, FOLLOWED 
BY THE RANGE FROM THE UNIT AT WHICH THE TARGET WAS 
DES TROYED 


"PENETRATION" IS PRINTED WHEN A CRUISE MISSILE APPROACHED 
WITHIN & KMs OF THE BATTERY. IN THIS INSTANCE THE 
UNIT IS CONSIDERED PENETRATED AND DESTROYED 


Table 3=-3.2 
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IV. DATA ANALYSIS 


Pee © URPOSE 

The purpose of this chapter is to demonstrate how the TI 59 may be 
used to analyze data. Rather than use assumed or contrived data, an 
actual experiment was conducted by the authors for illustrative purposes. 
Strict requirements for random sampling were not met in conducting the 
experiment but, again, the purpose of this thesis is to demonstrate the 
capabilities of the TI 59 rather than to make inferences or draw broad 
conclusions from the experimental data. 

Before discusSing methods for data analysis, the scope, design, and 
methodology of the experiment will be presented in sufficient detail to 
make the data analysis meaningful. Presentation of the experiment will 
be followed by a detailed discussion of a TI 59 program designed to 
compute measures of central tendency and spread for sample data. Chapter V 
discusses a TI 59 program which may be used to make statistical in- 


ferences using the same experimental data. 


Pee linESS EXPERIMENT 
ime scope Of the Experiment 
The experiment was conducted to test the hypothesis that 
different physical conditioning programs result in different levels of 
physical fitness. Six different conditioning programs were evaluated 
using five tests. The scope of the experiment was limited to testing 
the strength and endurance of selected upper body muscles, together with 


Overall cardiovascular fitness. A completely comprehensive fitness 


od 





evaluation would also include lower body strength and endurance as well 
as muscular flexibility and agility. Other factors such as diet, use of 
alcohol, tobacco and other drugs would also be requisite. This experi- 
ment was limited to the examination of thirty-six male subjects by the 
two authors to determine cardiovascular efficiency, bicep strength, 

bicep endurance, pectoral strength and pectoral endurance. Thirty of 

the subjects were military officers attending the Naval Postgraduate 
School. The other six subjects were weightlifters who trained at 
Bailey's Gym in Seaside, ae The subjects varied in age from 
nineteen to thirty-seven but were predominantly in their early thirties. 
While it is recognized that strict requirements for random sampling 
requisite for statistical analysis were not met in conducting the experi- 
ment, it should be pointed out that there is every reason to believe 

that the subjects examined were representative of the population from 
which students are continually drawn for the institution. Strict inference 
to any specified population will not be made, but as mentioned earlier 
the purpose of this discussion is to illustrate the use of programs de- 
veloped for the TI 59, not to make inferences from the data gathered in 
this particular experiment. 

Since the terms, strength, endurance, and cardiovascular 
efficiency are subject to a myriad of interpretations, the definitions 
used for this analysis are presented before proceeding to discuss the 
test methodology. Muscular strength is defined as the ability to 
exert maximum force against an object, while muscular endurance relates 
to the ability to exert force which. is not necessarily maximal over an 


extended period of time. Cardiovascular efficiency relates to how 
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well the heart, lungs and blood vessels work together without 
strain [ Ref. a : 

Before analyzing the results a thorough explanation of the 
experimental design and methodology is provided. 

2. Experimental Design 

The subjects were divided into six categories (each consisting 
of six individuals) based on their exercise programs. The following six 
categories were used: Category I - Individuals who had participated in 
no exercise over the last six months; Category II - Individuals who ran 
more than thirty miles per week and who did so for at least the last 
Six months; Category III - Individuals whose exercise consisted solely 
of lifting weights, whether for power or building physique, and who 
lifted at least twice a week; Category IV ~ Individuals who lifted 
weights at least twice a week, who ran twenty or more miles per week 
and who did both the lifting and running for at least the last six 
months; Category V - Individuals who lifted weights at least twice a 
week, who ran between ten and nineteen miles per week, and who did 
the lifting and the running for at least the last six months; and 
Category VI - Individuals who did not run or lift weights but who 
participated in basketball, racquetball, bicycling, or judo on a regular 
basis. 

Five tests were administered to each subject in each category. 
Cardiovascular efficiency was measured using the Pipes Test for Cardio- 
vascular Health, which consists of the following seven steps: 

1. Have the subject lie quietly on a pad for ten minutes. 


2. Take a pulse reading for the resting heart rate. 


a9 





3. Have the subject sit on a chair with his arms folded 
across his chest. 


4, With his arms folded, have the subject stand up and 
sit down twice every five seconds for three minutes. 


5. Take a reading for the heart rate immediately after 
sitting down at the conclusion of the three minute 
exercise. 

6. Take readings at 30-second, 60-second, and 120- 
second intervals after the exercise to measure 
recovery heart rate. 

7. Using a table developed by Pipes and the five heart 
rate readings determine the cardiovascular health 
score. 

This test, developed by an exercise physiologist, j Ref. a] is 
based on the premise that the heart rate describes an individual's 
fitness in three areas: how much. oxygen he needs, how much blood his 
heart has to pump to meet the oxygen need, and how hard the heart works. 
Individuals with a low level of fitness do not extract oxygen efficiently 
from the blood so must pump more blood, resulting in a higher heart rate. 
In conducting the test, the same examiner took the pulse readings at 
the different intervals using the radial artery in the subject's wrist. 
The heart rate was monitored for ten seconds then multiplied by six to 
obtain the number of heart beats per minute. 

The muscular strength test used for biceps was the maximum 
standing curl the individual could accomplish using a barbell. The 
proper technique was demonstrated to each subject by the examiner and 
lifts where the individual "cheated" by swinging the weight or arching 
his back were not counted. Bicep endurance was measured by the number 


of curl repetitions the subject performed with a 55 pound weight. These 


two tests were predicated on the generally accepted basis that curling 
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is the primary bicep exercise and that low repetition, high resistance 
exercises are best for developing endurance [ Ree. 3] | PLE tCy=f£ive 
pounds of weight were selected prior to the experiment a5 a low 
resistance weight well below each subject's strength capability and 
therefore in compliance with. expert opinion that the force used for 
endurance testing should be considerably below the individual's static 
force capability Ref. 7| : 

In a similar fashion, the bench press was used to test pectoral 
strength and endurance. Pectoral strength was measured by the maximum 
weight that the individual could bench press, while pectoral endurance 
was measured by the number of bench press repetitions he could perform. 

Each subject was tested using an identical sequence of events. 
Initially each person was given an instruction sheet (table B-1l in 
Appendix B) which explained the purpose of the experiment and defined 
those attributes to be measured, i.e. muscular endurance, muscular 
strength, and cardiovascular health. The subject was then asked to 
complete a questionnaire (table B-2, Appendix B) concerning certain 
aspects of his medical history. Each man was instructed to stop the 
testing if he felt any significant level of pain. He was then asked 
to complete a form disclosing his name, age, weight, and height 
(table B-3, Appendix B). Next, one of the examiners questioned the 
subject concerning his exercise program over the last six months and 
made a subjective judgement as to which of the six categories he 
belonged in. 

After this administrative procedure was completed, the actual 


testing was begun with the bench press test described above. (The 
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Same examiner tested each of the thirty-six subjects in both the strength 
and endurance exercises in order to minimize any variance due to test 
administration.) The subject was shown how to do the bench press, 
allowed to practice once if desired, and then tested for the maximum 
number of repetitions he could perform with 100 pounds. The number 

of repetitions was recorded and the subject was allowed a three minute 
rest before being tested for his maximum bench press. The examiner 
estimated the amount that each subject could bench press and set up 

the weights accordingly. All adjusting of weight was done by the 
examiner so that the subject's lift capability was not degraded. The 
man was then asked to bench press the weight set up for him. If he 

was able to make the lift ten pounds were added and he was asked to 

try again. If he failed the second attempt he was given credit for 
five pounds less than he attempted. For example, if an initial attempt 
of 165 pounds was successful and a subsequent attempt of 175 pounds was 
missed, then the score was recorded as 170 pounds. 

After another three minute rest the subject was tested on the 
number of times he could curl 55 pounds. This was followed by another 
three minute rest before testing for his maximum curl capability. Once 
again, the man's maximum lift was estimated by the examiner and all 
adjusting of the weight was done by the examiner. 

Following the four lift tests, the Pipes Cardiovascular Health 
Test | Ree. a | was administered by the second examiner and the subject's 
testing was completed. Before examining the test results, a procedure 
for scoring the tests was requisite. Accordingly, the scoring procedure 


explained in the following discussion was decided upon. 
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3. Scoring Methodology 


The heart rates recorded during the Pipes Cardiovascular Test 
were scored using table B-4 in Appendix B. Each subject's score for 
resting heart rate, heart rate immediately after the exercise and heart 
rate at the 30-second, 60-second, and 120-second intervals was aggregated 
to a total score ranging from zero to one hundred. This score was then 
used as the measure of cardiovascular fitness for comparative analysis. 

In order to compare muscle strength among the subjects and 
among the categories, it was necessary to adjust each subject's lift for 
varying sizes and body structures. Accordingly, each man's maximum curl 
and maximum bench press were divided by his body weight, resulting in 
an adjusted score for each lift. These two adjusted scores were then 
added together to yield an upper body strength measure. For example, 
let Sj3 be the strength measure for the third subject in Category I 
(where the first subscript indicates the category and the second indicates 
the subject within the category). The following formula may then be used 


to obtain the strength score for the third subject tested in Category I: 


= maximum bench + maximum curl 
body weight 


As discussed earlier, muscle endurance was measured for the 
Same two areas tested for strength - the biceps and the pectorals. The 
bicep endurance was measured by the maximum number of curl repetitions 
performed with 55 pounds, while the pectoral endurance was measured by 


the maximum number of bench presses accomplished with 100 pounds. As 
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in the case of strength, an adjustment was made for the subject's body 
Wemeanct. In the case of endurance, however, the amount of weight lifted 
(which was 55 pounds for the curl and 100 pounds for the bench press) 
was divided by the subject's body weight and then multiplied by the 
corresponding number of repetitions lifted. These two scores were then 
summed as the endurance index. For example, let E13 be the endurance 
score for the third subject tested in Category I. The following formula 


then obtains: 


= _ 100 ibs Ss * (number of bench press 
13 body weight of subject repetitions) 
Es 25 lbs * (number of curl 
body weight of subject repetitions) 


4. Test Results 
Tables B-5 through B-10 in Appendix B reflect the results of the 

experiment for each of the six categories tested. For example, table B-=-5 
depicts the age, weight, cardiovascular score, adjusted strength score 
and adjusted endurance score for each of the six subjects tested in 
Category I. Appropriate references are made at table B-5 for the develop- 
ment of the final cardiovascular, strength and endurance scores. The 
scoring methodology section of the chapter provides a detailed explanation 
of the rationale and methodology for deriving these scores. 
C. TI 59 PROGRAM FOR DATA ANALYSIS 

Having developed the experimental design, the scoring methodology, 


and the test results it is now possible to analyze the data. Measures of 
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central tendency and spread will be used to illustrate an application of 
the TI 59 in analyzing sample data. The measures of central tendency 
used for this illustration are the mean, X¥ , and the median. The 
measures of spread used are: the sample variance, s¢; the standard 
deviation, s ; the mean absolute deviation, MAD; the mean squared 
deviation, MSD; the root mean squared deviation, RMSD; and the range. 
A TI 59 program will now be described in detail which computes these 
measures followed by an example applying the program to the results of 
the fitness experiment. 
ie Li 59 Capabilities 

The TI 59 has been hard-wired to calculate the sample mean and 
variance as well as MSD. As described in the TI 59 Personal Programming 
Manual [Ref. 12] if each datum is entered into the calculator followed 
by pressing the > + key, the calculator will sum each data entry, Xj, 
into register one, sum the squares of x; into register two, store the 
number of data entries in register three, and calculate xX, s?¢, and MSD. 
fayodeftinition, s? = 1l/n-1 >) (x1 - x) ¢ is the unbiased estimator for 
g* while MSD = (n-1)s¢/n is the maximum likelihood estimator.) 
Pressing the xX key will yield the mean, INV X will display s? and 2nd Op 
ll will display MSD. If these are the only measures desired then 
utilization of the ».+ key is the most expedient method of obtaining 
them. The TI 59 statistics module has a program (Program 03) which 
computes these same measures as well as the middle value (MIDVAL). 
Additionally, Program 03 stores each data entry beginning with register 


31. Program 03 also computes a number of other quantities not germane 
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to an analysis of the data gathered in the experiment discussed pre- 
viously in this chapter. Since this results in a slightly longer run 
time for each computation, a program has been written by the authors 
which exploits the hard-wire capabilities of the > + key, computes MAD 
and range in addition to the other measures discussed, and stores the 
data for recall or transformation if desired. In addition, this pro- 
gram may be used with the TI 59 Master Module if the Statistics Module 
1s not available. The following section describes the program in detail. 
2. Univariate Data Program 

In order to facilitate the description of this program a flow- 
chart (figure 4-1), has been included at the end of the chapter. 
Comments in the paper are keyed to figure 4-1 by numbered circles for 
easy reference. 

The program is initialized by pressing E' (figure 4-1,G)). 
Initialization entails clearing all of the data registers, lowering 
flag 0 the purpose of which will be addressed later, and storing 31 in 
register 30. Register 30 is used as a post office for indirect address- 
ing. In this particular program this means that data are stored in the 
register indicated by register 30. For example, after initialization, 
Begrster 30 contains 31. The sequence x), STO 2nd IND 30 will result in 
xy being stored in register 31. When the initialization routine, 2nd E&’, 
is complete the display will contain the value 31. Each datum mav now be 
entered successively followed by pressing A. The routine at Label A 
begins by storing x; in a working, register 18 (figure 4-1,(2)). 


x; 1s then stored permanently beginning in register 3l. xy 1s stored 
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in register 31, then register 30 the indirect storage address is incre- 
mented by 1 so that X5 wil ebemeeercd ineregisters32, x. in register 33 


and x, in register 31+n- 1. A total of sixty-nine entries may be 


n 


made uSing registers 31 through 99 for data storage. Registers 0 
through 29 are used to make the requisite computations of central 


tendency and spread. After each datum, x: 


Ane has been stored, the program 


checks to see if flag 0 is raised (figure 4-1,(3)). Tete lagacs Ls 
raised this indicates that a data entry has been made previously, ie., 


Ege cirrent x; is not xX) - In this event the program skins to Label x. 


If flag 0 is not raised, ie., the current x: is x 


4 17 then x, is recalled 


from the working register, register 18, and stored in register 12 as 


the minimum x. and register 13 as the maximum x Future entries may 


i! 1 
then be checked against register 12 to determine which value is lower. 


If a current x; 1s lower than the value in register 12 then it will 


1 
replace it as Xnin: Similarly, subsequent entries may be checked 
against register 13 in order to retain Xnax- Peter Storing x) in 
register 12 and in register 13 the program internallv calls the key 
which, as disscussed previously, will sum Xs into register Ol, sum 
a into register 02, sum the number of entries into register 03 and 


compute xX, s¢, and MSD. Flag 0 is then raised so that subsequent 


entries will skip to Label x and replace or as appropriate. 


max 
The program then recalls the number of entries, n, into the display and 
stops awaiting the next entry (figure 4-1,(4)). The second entry x5 will 
now be stored temporarily in the working register, register 18, and 


permanently in register 32. The indirect addressing register, register 


30, is incremented by 1 for the next entry and the program then checks 
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to see if flag 0 is raised. Since this is not the first entry, the flag 
will be raised causing the program to skip to Label x (figure 4-1, (s) re 
The first step under Label x is to recall nin from register 12 and 
store it in the test register Rp» For this particular iteration, x 
will be in register 12 since the first entry was both the maximum and 
the minimum value processed as described earlier. The program then 
recalls the current x; (x, in this instance) from the working register, 
register 18. The display value, X54 is checked against the Rp value, 
Xs fo see if the display value is less than the Rp» value. If so the 
program skips to step 57 where Xo is stored in register 12, replacing 

x, as the lowest data entry (figure 4-1,(6)). If the display value is 
not less than the Re value then the program recalls Xnax from R13 and 
stores it in Rp: Xj is recalled from the working register, register 18, 
into the display. This time the program checks to see if the display 
Werme, xX, 1S greater than the Rp value, eae If so, the program skips 
to step 62 where x; ils stored in register 13 as the new ae (fig. 4-l, 
@). The program then computes the MIDVAL by recalling Sa from 
register 12 and Xmas from register 13, summing them and dividing by 

2, and storing in register 14. Next the range is computed by sub- 
tracting a from x nax’ The range value is stored in register 15 
(figure 4-1,(8)). The program then loops back to the + key to 

compute the mean, variance and MSD, (figure 4-1,(9)). The number of 
entries is recalled and displayed awaiting the next entry. This process 


1s repeated until all of the data have been processed. 


The outputs of the program may be recalled as shown in table 3-ll 
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in Appendix B. The mean is displayed by pressing x, the variance by 
pressing Inv X, and MSD by pressing 2nd op 11. The lowest data point, 
Xmin’ May be discovered by recalling register 12 while the highest data 
entry, Xnax’ May be recalled from register 13. Recalling register 14 
will display the MIDVAL and the range may be found by recalling 
register 15. Each of the original data entries may be recalled if 
desired beginning with x, in register 31. MAD is computed by pressing 
2nd A' which calls a different subroutine. This subroutine recalls x 
which was computed under Label A and stores it in Register 16. The 
number of entries, n, is recalled from register 3 and stored in 
register 7 to be used as a decrement register. Register 20 contains 
31 which is used to indirectly address the datum which have been stored 
beginning with register 31. The program recalls each x5 using register 
20 and subtracts x. The absolute value of the difference is summed 
anto register 19. The program does this successively for each x; until 
the decrement register, register 20, is equal to zero indicating that 
each Xs has been processed. The sum of the absolute values of the 
deviations from the mean is recalled from register 19 and divided by 
n which is recalled from register 3. This value, the mean absolute 
deviation is displayed completing the subroutine A' processing. All 
the values discussed earlier are still intact and may be recalled if 
needed. Table B-12, Appendix B, is a program listing for the uni- 
variate program. 
D. APPLICATION OF THE TI 59 UNIVARIATE PROGRAM 

The cardiovascular scores for Category I provide a ready example for 


the use of the univariate program to calculate measures of central 
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tendency and spread. After the program card has been read in, the 
program is initialized by pressing 2nd E*. The cardiovascular scores 
for Category I (table B-5) are entered into the calculator as follows: 
memeeenA; 5S, A; 58, A; 44.5, A; 33.5, A; 40, A. The instructions 
contained in table B-ll may then be used to obtain the desired statistics. 
For this particular example: 2nd x yields the mean, 48.4; RCL 14 dis- 
plays the MIDVAL, 45.75; RCL 15 displays the range, 24.5; 2nd A' yields 
the mean absolute deviation, 9; 2nd Op ll displays the mean squared 
Memration, 92.9; 2nd Op l1l\f/x , calculates the root mean squared 
deviation, 9.6; INV 2nd x recalls the standard deviation, 10.56; INV 

end xX, x* calculates the variance 111.5; RCL 12 displays the lowest data 
entry, 33.5 and RCL 13 displays the highest data entry, 58. To 
calculate the sample statistics for another category or for a different 
test the user need only push 2nd E' to re-initialize and then enter the 
relevant data. Statistics have been calculated for the age, weight, 
cardiovascular scores, endurance scores, and strength scores for each 

of the six categories. Tables 4-1 through 4-6 display these statistics. 
Rather than discuss each of these tables in depth, one example is pro- 
vided relative to the interpretation of the sample statistics. 

The cardiovascular mean for Category I, 48.4, indicates average 
cardiovascular fitness using the Pipe's test which is based on a scale 
from 0 to 100. Three measures of spread (standard deviation 10.56, 
root mean squared deviation, 9.6, and mean absolute deviation, 9) are 
approximately equal to ten, a rather high variability in this case. 


The range, 24.5, also indicates that the data are quite spread out. 
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Xin of33.5 and See of 58, the bounds of the sample data indicate that 
the cardiovascular fitness of sedentary people varies from poor to 
average. 

Inferences, subject to the sampling limitations already discussed, 
may also be made about the strength or endurance of sedentary people 
using the data from table B-5. Similarly, the statistics for the. 
other categories may be used to make inferences about the strength, 
endurance, or cardiovascular fitness of those who run over thirty miles 
per week (Category II) or those who lift weights (Category III) or 
any of the other three categories. Programs have also been written wor 
the TI 59 which allow a user to develop confidence intervals for these 
sample statistics Ref. 14| : 

The next chapter will discuss a program for one factor analysis 


of variance and then apply the program to the fitness data to illustrate 


Statistical inference with the TI 59. 
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Vo PAT USTY CAL INP ERENCE 


A. PURPOSE 

The purpose of this chapter is to illustrate a method of statistical 
inference using the TI 59. As in chapter IV the intent is to demonstrate 
the capabilities of the TI 59 rather than to emphasize statistical prin- 
ciples. One-way analysis of variance (ANOVA) will be used to demonstrate 
statistical inference using the fitness data discussed in the preceeding 
chapter. The variations of the underlying populations represented by the 


Six categories are assumed to be unknown but equal for this illustration. 


B. TI 59 PROGRAM FOR ANALYSIS OF VARIANCE 
In testing the hypothesis that the population means for each of the 


Six test categories are equal, H, is typically rejected if the F ratio 


0 
exceeds the critical F value in a standard table for the desired test 
level (typically 5%). Alternatively, using the TI 59 to its full 


advantage, prob-value may be used to test H Prob-value is a method of 


0° 
testing whether or not the null hypothesis is supported by the data. In 
the case of the F ratio, prob-value is the probability that the F ratio 
would be as large or larger than the value actually observed if Hy were 
true. This is the right hand tail area, Q (f), where 

OC yes Pyrat > o£) 
Prob-value has the advantage that analysis is not restricted to arbi- 
trarily established test levels such as 5% or 10% or to use of standard 


published tables. The TI 59 Statistics Module has an F distribution 


program (Program 22) which computes the tail area of an F curve where 
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the curve is defined by the degrees of freedom in the numerator and the 
denominator. A series expansion is used to approximate the integral to 


determine Q(f) [Ref. 10| - tte ais true, indicating that all of the 


0 
observations are from the same normal population, then the prob-value, 


O(£), will be large. Conversely, if H, is false then the prob-value 


0 
will be small. If the prob-value is sufficiently small (as determined 
by the decision-maker) then Hy is rejected and the conclusion is formed 
that there must be a difference in the population means somewhere 
[Reé. 14] | 
While a classical test or prob-value may facilitate rejection or 
acceptance of the null hypothesis, no insight is provided as to which 


means differ, given that H. is rejected. There is an efficient method 


0 
developed by Sheffe [Res. 14 for computing confidence intervals for 
the difference between means. If the physical fitness example discussed 


earlier is used, then Sheffe's development may be used to make the 


following statements with 95% confidence: 


VRE Svar 


(4, — H,) V (r- Es 5 n, : n; 
H-4) = (K-%) + Vr-RSVE, + 4, 


{| 
“oN 
<| 
<| 

Ww 
ee” 


(H, — H,) 


| 
>! 
us 
+ 


Where F 05 is the critical value of F which leaves 5% in the upper tail, 


Sp is the square root of the pooled variance, r is the number of means 
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compared, andn is the size of each of the samples. In the fitness 
example this equates to making confidence statements about the difference 
in fitness between the six categories where r is six and n is 6 for each 
of the samples. To facilitate multiple comparisons a contrast of means 


is used. This contrast may be written as: 


ye Cc U. where ye C. =) 


It is then possible to develop the following formula which includes all 


possible contrasts with 95% confidence: 


Uoa= Lox *Vimje, Sp VE 


If the value O is included in a confidence interval then there is no 
basis for believing that the population means differ while if O is not 
included then the conclusion is drawn that the means do in fact differ. 

A program has been developed by Dr. P.W. Zehna Res. 14 for the 
TI 59 which computes the elements of an ANOVA table (table B-15)to 
include the F ratio and prob-value discussed above. Basically the 
program exploits the TI 59 F distribution program for determining 
prob-value (Program 22) after using Program 15 of the Statistics Module 

Res. 10] to calculate the F statistic. The program then uses Scheffe's 

multiple contrasts to determine which population means differ given 
that the null hypothesis is rejected. 

A flowchart (figure 5-1), user instructions (table B-13), and a listing 
of the actual program steps (table B-14) are provided to facilitate 
description of the ANOVA program. The program takes data input by rows 


and outputs the elements of an ANOVA table (table B-15) sequentially 
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as indicated by the number in each block of the table. A row of data 
constitutes a sample as in the example at table B-16. 

The program begins by using Program 06 of the Statistics Module to 
enter the data. After initialization with 2nd E' each 8 is entered 
followed by pressing Label A. When one complete row has been entered a 
press of 2nd B‘' causes the calculator to compute the row or sample x and 
a press of 2nd C' results in computation of the MSD. These two steps 
must be performed after each row has been input so that the calculator 
will know when a new row is being entered. When all of the data have 
been entered using this scheme (table B-13) the sequence RST, A begins 
the ANOVA Table calculations. The first step under Label A is to call 
Program 15 of the TI 59 Statistics module which computes the F ratio. 
In the process of computing the F ratio the other elements of the ANOVA 
table (table B-15) are computed and stored except for the prob-value. 
To fill in the values for the ANOVA table all that is required is 
successive pushes of R/S as indicated in table B-13. For example, the 
first R/S displays the degrees of freedom for the numerator while the 
fourth R/S displays the degrees of freedom for the denominator. The 
program essentially recalls and displays the calculations of Program 15 
of the Statistics Module to build the ANOVA table. To compute prob- 
value the program internally calls Program 22 of the Statistics Module. 
The user need only press R/S as indicated in table B-13 which causes the 
calculator to recall the degrees of freedom for the numerator and the 
denominator used in Pgm 15 and transfer them to Program 22 to define 


the F Distribution. The F statistic calculated in Program 15 is then 
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recalled and transfered to Program 22 resulting in O(f), the probability 
that F> £. This prob-value may then be used to accept or reject the 
null hypothesis. If the null hypothesis is accepted then the analysis 
is completed. However, if Ho is rejected, the next step entails the use 
of Scheffe's contrasts to determine which means differ. 

To use the ANOVA program (table B-14) for posterior contrasts with 
Scheffe's formulas the user initializes the routine by pressing 2nd E'. 
Then cj, Xs and n, are entered for each row as shown in table B-13. The 


eS are the coefficients used to determine which means are contrasted 


as discussed previously. To contrast u, and u5, c, = l, cy = -l and all 

other c; = 0. To contrast uy and Us, Coa Py Go eet ended tones 

c,; = 0. As these data are input, the program uses a ‘'loop' to calculate 
& 


(=) which is stored in register 03 and >. ¢ xX: which is stored 
\ 
in register 06, (figure 5-1). Register 04 is used as a counter to 


display the number of row entries. After each ¢;, 


Xs and ns entry, 
register 04 is incremented by one. The program then transfers to 

Label x, displays the running count of row entries and stops pending the 
next entry. After every cj, x; and n, have been processed, a Eritical 
value of F with degrees of freedom r-l and n(r-1) is entered for the 
desired test level followed by 2nd A'. The program recalls the degrees 
of freedom for the numerator, r-l, from register 14 and the pooled 
varlance S. from register 29. The product (r-1l) F2)S." is formed and 
multiplied by the contents of register 03 ’. (=) . The square root 


of this product is stored in register 05. This value is then added to and 


subtracted from the contents of register 06, Dic. <5 to form the desired 
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confidence interval. The lower bound is displayed after the use of 
2nd A' and the upper bound may be recalled by pressing R/S (table B-13). 
An example will now be provided uSing this program to test for 


differences between population means for the fitness experiments. 


C. APPLICATION OF THE TI 59 ANOVA PROGRAM 


The null hypothesis, H., may be stated as - there is no difference 


Oo! 
in the cardiovascular fitness of those who do no exercise (Category I), 
those who run in excess of thirty miles per week (Category II), those 
whose exercise consists solely of lifting weights (Category III), those 
who lift weights and run in excess of twenty miles per week (Category IV), 
those who lift weights and run between ten and nineteen miles per week 
(Category V), and those who do not run or lift weights but participate in 
other activities such as basketball, racquetball or bicycling (Category 
VI). Table B-35 reflects the cardiovascular score for each of the thirty- 
six subjects tested by category as well as the mean for each category. 
The null hypothesis that there is no difference between these category 
means may be tested using the ANOVA program with each of the categories 
Senetttuting a row for input. 

After the program card has been read in, the ANOVA program is used 
by first calling program 06 of the statistics module to enter the 
data. After initialization with 2nd E', the data for each row are 
entered followed by A. For the cardiovascular scores (table B-35) the 
first row, Category I, would be entered as follows: 56.5, A; 58, A; 
44.5, A; 33.5, A; 40, A. Once the row data have been entered 2nd B' 


is pressed to display the row mean, 48.4, followed by C' which displays 
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the row MSD, 92.9. (The row mean must be recorded for use in posterior 
contrasts). The data are then entered in a similar fashion for the 
remaining five rows (Categories II through VI). Once all of the data 
have been entered, RST is pressed to return the calculator pointer to 
the ANOVA program. A is then pressed resulting in calculation of the 
ANOVA Table entries. The ANOVA entries are recalled with sequential 
presses of R/S. Table B-13 discussed earlier contains detailed 
instructions on the use of the ANOVA program. Table 5-1 depicts the 
ANOVA calculations for the cardiovascular scores of the six fitness 
test categories. The prob-value of .00027 is sufficiently small to 
cast doubt upon the null hypothesis that there is no difference in the 
cardiovascular fitness among the six categories tested. 

As discussed previously, the prob-value tells how credible the null 
hypothesis is but it does not tell which categories differ given that 
Eaere 1S cause to reject Hg. However, confidence intervals may be 
established for contrasts between the categories using Scheffe's 
formula which is programmed in the ANOVA routine Ref. 13| 2 oy ene 
cardiovascular example, the cardiovascular fitness of the sedentary 
subjects (Category I) may be contrasted to the cardiovascular fitness 
of the runner (Category II) as a demonstration of the program. 2nd E' 
is pressed to initialize the contrast routine followed by cj, x, and 


mn; for each of the two rows. For Categories I and II the entries are: 


iL R/S 48.4 Ryo 6 R7S 


-1 R/S 83.4 R/S 6 R/S 
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The appropriate F percentile is entered followed by A' to generate 
the desired confidence interval. To display a 95% confidence interval 
for the difference in cardiovascular fitness between Categories I and 
il an F Bee centile of 2.53 (where there are five degrees of freedom in 
the numerator and thirty degrees of freedom in the denominator) is used 
resulting in an interval from -65.5 to -41.4. Since 0 is not included 
in the interval it is reasonable to conclude that there is a difference 
in the cardiovascular fitness of the two categories. Table 5-2 con- 
tains the results of contrasting each of the six fitness categories. 
Only four contrasts result in the conclusion that there is a difference 
between the categories with 95% confidence: Category I - Category II 
(-65.5, -41.4); Category I - Category IV (-65.5, -4.4); Category II - 
Category III (.6, 61.8) and Category III - Category IV (-61.7, -6.Q). 
These results indicate with 95% confidence that there is a difference 
in the cardiovascular fitness of those who run more than twenty miles 
per week (Categories II and IV) and those who do no running at all 
(Categories I and III), at least for those subjects examined. 

The ANOVA program has also been applied to the strength and 
cardiovascular scores resulting from the experiment. Table B-36, 
Appendix B, reflects the strength scores of each of the thirty-six 
subjects by category. The ANOVA results are contained in Table 5-3. 
The prob-value of .12 x 107° indicates that the null hypothesis that 
there is no difference in the strength of members of the different 
categories should be rejected. Further analysis with posterior con- 
trasts is necessary to see which categories differ. Table 5-4 contains 


the results of posterior contrasts with an F percentile of 2.53 for 25% 
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confidence with five and thirty degrees of freedom for the numerator and 
denominator respectively. Unlike the cardiovascular contrasts there 
are a number of differences in the strength results. Categories III 
and IV, which were composed of the most ardent weightlifters, differs 
from Categories I and III but not VI. These results are not surprising 
in that they confirm the hypothesis. that different training programs 
result in different levels of fitness. In this instance where fitness 
is defined as strength, those who trained for strength were in fact 
stronger than those who did not. Again, without attempting inference 
to a larger population, these results may be used to gain insight into 
the probable differences that might be tested in a more appropriately 
designed experiment. 

Table B-37, Appendix B depicts the endurance scores for the thirty- 
Six subjects by category. Table 5-5 reflects the results of using the 
ANOVA program with the endurance scores as input. Once again, the 
prob-value of .137 x 10-© indicates rejection of the null hypothesis 
(in this case that there is no difference in the upper body endurance 
of the members of the six different categories). The posterior con- 
trasts (table 5-6) indicate that the weightlifters (Categories III, 
IV, V) differ from the non-weightlifters (Categories I, II, VI) in 
Mpper body endurance with 95% confidence. This also supports the 
hypothesis that different training programs result in different levels 
of fitness, subject again to the sampling restrictions previously 


discussed. 


oa 





D. SUMMARY 

While the results of the fitness experiment are interesting, the 
purpose of this analysis has been to demonstrate a statistical applica- 
tion of the TI 59 and not draw inference to a hitherto undefined popula- 
tion. The univariate program was used to calculate measures of central 
tendency and spread for the Category I Cardiovascular scores. The 
ANOVA program was used to test for differences in strength, endurance, 
and cardiovascular fitness among the six test categories. In both 
instances meaningful but guarded inferences were drawn from the test 
data. 

The capabilities of the TI 59 in real world statistical analysis 
are impressive. The analyst can conduct sophisticated analysis of 
good-sized samples unconstrained by access to large computers. Using 
programs such as those demonstrated in this chapter the analyst need 
not even learn a programming language. All that is required to compute 
an F ratio or prob-value, for example, is the ability to follow simple 
users' instructions. While there are certainly samples whose size 
preclude the use of the TI 59, there are a pletheora of samples which 
can be analyzed more conveniently and just as efficiently at home or 


at the office using the TI 59. 
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APPENDIX A 


COMMENTS ON SELECTED LABELS 


LABEL COMMENT 
A Directs action to compute glimpse probability of detection 


for both ICWAR and IPAR 


[| Computes glimpse probability of detection for ICWAR only 
based on target range 


| 


Computes glimpse probability of detection for IPAR only 
based on target range 


E and Lnx Determines if target is less than 40 KM from unit and 
directs processing accordingly. 


Dp: Calculates target range after elapse of target lock-on time. 
fopR) EXC and INT Increments range of IHAWK missile after firing 

FIX Prints 24 (no firing capability) 

DSZ Prints 14 (cruise missile identified as friend) 

NOP Prints 66 (IFF is nonoperational) 

OP Prints 1 (target detected) at range (KM), begins engagement 

sequence 

(SBR) RAD and y* Increments simulated air battle time clock (Rj). 
(SBR) GRAD and x Computes cruise missile rate of approach and incre- 


ments target range (Roo) - 
(SBR) ENG Random number generation - normal distribution 


(SBR) P—»R Random number generator - uniform distribution 


PRD Prints 23 (no detection capability) 
SIN Prints "KILL" and range of target kill 
LOG Prints "PENETRATION" 
(SBR) WRITE AND (‘TT Prints "PENETRATION" if target is 8 KM or less 


mremelinial te. 


Le 





START 





ag 
ePAR OP? 


NO 


SET FLAG 3 
PRINT 3 


LBL DEG 


“SET FLAG 4 
“PRINT 4 
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LBL PRD 
BATTERY 
NONOPERATIONAL 
NO DET CAP 


NO 


€os 
IS 


TPAR OP: 


YES 


"SET FLAG 5 
*PRIND ge 


EEE CG 
TS 
TARGET DET? 


YES 





NO TARGET es 












DETECTION 
LBL LOG 
PENETRATION 
NO 
LBL RAD 
SLOr 
CLOCK 
LBL GRAD 
INCREMENT 
TARGET RANGE 
YES 
TARGET VED 
GTO O00 RIEND? 
NO 
LBL A 





YES 
ICWAR + IPAR (a) 


PommAG 53 ob. ¢ 


NO 


Is LBL_B 


YES 
ONLY ICWAR 
IS FLAG 4 SET? OP 


NO 


LBL — 


1LO4 










Sei Or 
SET FLAG 1 
PREN = 


RCL TARGET RANGE - PRINT 







YES LBL NOP 


SET FLAG 6 
PRINT 6 













e000. Fas 2 


IS 
TARGET 
FRIEND? 











Gio LEL STO 








YES 


PRINT 10 


BL 3, LBL ADV 















NO NO 
A SECTION B SECTION 
PRINT 18 OP? 
YES 
PRINT 17 LBL FIx 
vas BATTERY 
NONOPERATIONZ 
NO TRK CAP 
B SECTION 
STOP 
BL LIST 
ASSIGN ASSIGN ASSIGN 
ee TO A TO B 


SET S 
SET FLAG 0 gPRINT 9 FLAG 7 7 PRINT 7 AG 8 Y PRINT 8 


CALCULATE 





CALCULATE 
EOCK=ON TIME 





CALCULATE 
LOCK-ON TIME 


SBR ENG 









LOCK-ON TIME 





NOS 








ES FLAG O° SET 


LBL D' Lane a. 


CALCULATE TGT 
MOVEMENT DURING 


CALCULATE TGT 
MOVEMENT DURING 
EOCK-ON TIME= 

BOTH SEC 


LOCK-ONe TIME= 
ONE SEC 
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Ig LBX xX 


TARGET 
£3 KM? 
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eS YES 
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PRINT 12 










INCREMENT 
CLOCK AND 
TGT RANGE 






Aes LON 
MISSILE AWAY 


B SECTION 












MISSILE AWAY 


INCREMENT 
CLOCK 
TGT RANGE 

> MSL RANGE 


LBL GRAD 
LBL RAD 
LBL INT 
LBL T 










LBL WRITE 
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YES 
ee TARGET 
See <8 KM? 
NO 


LBL SIN 
STOP YES 
STOP 
YES Ret CERIN Me Ler. SiN 
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YES 
LBL GRAD 
LBL RAD 
INCREMENT 
CLOCK AND 
TARGET RANGE 
LBL EE 
MISSILE AWAY PRINT 1l 
LBL E* 
LBL GRAD 
LBL RAD 
ELBE er 
LBL INT 
LBL 7 
=O 
NO 
LBL CLR 
NO YES 
SLOF 
YES 
LBL SIN 
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SIMULATION PROGRAM 
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APPENDIX B 
INSTRUCTION SHEET 


This is a physical fitness test. 
I. WHAT WE ARE EXAMINING AND WHY 

The following tests will look at the complex systems that make up 
the network of health and fitness, First we will test your upper body 
endurance and muscular strength. Next we will examine your cardiovas- 
cular health. This test is a component of a research project being 


conducted by graduate students of the Naval Postgraduate School. 


IT. MUSCULAR ENDURANCE AND STRENGTH 

Muscular endurance is often synonymously and incorrectly used in 
place of muscular strength. Muscular strength is the ability of your 
muscular system to exert maximum force against an odject or reSistance 
all at once, your ability to exert a maximum force a single time. 
Muscular endurance relates to the ability to exert force, not necessarily 
maximal, over an extended time period. As with all the components of 
fitness, these two concepts are interrelated but distinctly different 
from each other. Each concerns itself with particular capacities of 


fitness. 


ITT. CARDIOVASCULAR HEALTH 

Objective: To measure your heart-rate response to exercise. This 
test will determine the relative efficiency of your heart and circula- 
tory system. Your heart rate can be used to describe the fitness level 


of your body in three areas: how much oxygen you need, how much blood 
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your heart must pump to supply this need, and how hard your heart must 
work at this task. If your need for oxygen is not being fulfilled, your 
body is working in an inefficient manner. Consequently, more blood will 
have to be pumped through your circulatory system at a faster rate to 
get the oxygen to the muscles and organs that need it. The heart has 
the responsibility of satisfying your body's need for oxygen. It will 
have to beat more frequently to circulate the blood throughout your 
eyovem. tf your body works in an efficient manner, its need for 

oxygen is being fulfilled. 

Thus cardiovascular health relates to the ability of the heart, 
lungs, and blood vessels to work in unison without strain. Regardless 
of what the task is, whether physical or mental, the cardiovascular 
system should be able to handle it. When you have high levels of 
cardiovascular health you perform with more efficiency and you are 


more effective at what you do. 


TABLE B-1.2 


ZO 





QUESTIONNAIRE 


THIS IS A PHYSICAL FITNESS TEST 
READ THIS DOCUMENT CAREFULLY BEFORE THIS TEST BEGINS 


(1) This test is completely voluntary! You may decline testing now if 
you wish. You may stop at any time during this test and decline further 
testing. You are under no obligation to complete this test. 
(2) Before proceeding with this test, you should assure yourself and 
your tester that there have been no incidents in your medical history 
that would prohibit you from pursuing this testing. Your medical history 
is relevant to this physical fitness test. 
(3) Please answer the following statements: YES or NO 

(a) I have a heart related disease. 


(b) I have high blood pressure. 


(c) I often feel faint and suffer spells 
of dizziness. 


(d) I have recently or in the past felt pain, 
heaviness or pressure in my chest. 


(e) I have felt pain, heaviness or pressure 
in my chest when I walk uphill. 


(f) My doctor has advised me not to engage in 
physical exercise or physical activity. 


(4) I have read and fully understand this document. 


DATE SIGNATURE 
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DATA SHEET 
NAME 
AGE WEIGHT DATE 
PHYSICAL ACTIVITY CATEGORY: I II i ene + VL 
STRENGTH TEST: 
Maximum Bench Press 


Maximum Curl 


ENDURANCE TEST: 
Bench Press Repetitions With 100 Pounds 


Curl Repetitions With 55 Pounds 


CARDIOVASCULAR TEST: (Pulse Rate) 
Resting 


Immediate Post Exercise 
30 Seconds Post Exercise 
60 Seconds Post Exercise 


120 Seconds Post Exercise 
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CATEGORY I RESULTS 


SUBJECT NEE WEIGHT CARDIOVASCULAR STRENGTH ENDURANCE 
SCORE (1) SCORE (2) SCORE (3) 

fl 29 190 56.5 1.26 eet 

Z 33 155 58.0 ieee 1259 

3 36 170 58 .0 97 Dee 

4 31 205 Le. 5 1.19 270) 

» 37 7 poe ieate 10.1 

6 34, 160 40.0 1.16 11.0 


(1) See Table B-17 for derivation 
(2) See Table B-18 for derivation 
(3) See Table B-19 for derivation 
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CATEGORY If RESULTS 


SUBJECT AGE WEIGHT CARDIOVASCULAR STRENGTH ENDURANCE 
SCORE (1) SCORE (2) SCOnn 3) 

i 33 167 91.5 teeth 7 220) 

2 34 195 96.5 4 44 28.9 

3 35 180 55.0 dete 12e2 

4 By 160 93.5 1.44 1366 

5 28 127 e220 fe 6.5 

6 30 158 2220 1.414 6.6 


(1) See Table B-20 for derivation 
(2) See Table B-21 for derivation 
(3) See Table B-22 for derivation 
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CATEGORY III RESULTS 


SUBJECT AGE WEIGHT CARDIOVASCULAR STRENGTH ENDURANCE 
SCORE (1) SCORE (2) SCORE EGS) 

4 22 200 35.0 2 48 50.0 

2 23 202 56.0 25 62.5 

3 28 205 Ay. 5 2.49 Bae 

4 3, 220 86.5 22s 2 a7 

5 25 180 67.0 2.83 54.8 

6 29 225 24.5 2.00 oleae 


(1) See Table B-23 for derivation 
(2) See Table B-24 for derivation 
(3) See Table B-25 for derivation 
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CATEGORY IV RESULTS 


SUBJECT AGE WEICHT CARDIOVASCULAR STRENGTH ENDURANCE 
SCORE (1) SCORE (2) SCORE (3) 

1 28 200 96.5 230 47.5 

2 34 180 79.0 2.67 49.0 

3 32 165 85.0 2.00 43.7 

4 Bz 160 Fetal 2.19 Wee 

5 19 240 89.5 1.85 34.3 

6 26 180 Va 5 2739 Ly 6 


(1) See Table B-26 for derivation 
(2) See Table B-27 for derivation 
(3) See Table B-28 for derivation 
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CATECORY V RESULTS 


SUBJECT AGE WEIGHT CARDIOVASCULAR STRENGTH ENDURANCE 
SCORE (1) SCORE (2) SCORE (3) 

4 30 165 39.5 2.03 25.5 

2 24 ie 55.0 2.43 54.1 

© 19 ies 47.0 1.86 47.0 

4 34 79 GD 1.73 Sop 

5 35 205 59.5 des Zoe 

6 30 120 56.5 1.75 23.8 


(1) See Table B-29 for derivation 
(2) See Table B-30 for derivation 
(3) See Table B-31 for derivation 
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CATEGORY VI RESULTS 


SUBJECT ACE WEIGHT CARDIOVASCULAR STRENGTH ENDURANCE 
SCORE (1) SCORE (2) SCORE (3) 

a. ley 180 SI. fle 2a 

2 32 125 47 5 17556 i 

5 Bt 170 65.5 1.26 jee 

4 34 170 85.0 1.44 Pore 

5 33 175 86.5 1.34 11.4 

6 30 205 80.0 tS 17.9 


(1) See Table B-32 for derivation 
(2) See Table B-33 for derivation 
(3) See Table B-34 for derivation 
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Univariate User Instructions 





Initialize 
Enter data x. 
Repeat for each x5 
Recall Statistics one x 
INV 2nd & 
2nd OP ll 
RE be 
REE) 13 
RCL 14 
RCI, 
2nd 
RCL 
Recall data entered in 
Step 2 if desired. 
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ANOVA USER INSTRUCTIONS 





Select Program 06 2nd Pgm 06 
Initialize data base end b' 
Enter data for each row X43 A 
Reset pointer if more 
than 29 data entries 
are made D 
Gal culate ode el ON ond B' 
Calculate MSD for 
current row ene 
Revue tO ustep 9° tO 
enter next row data 
Return pointer to 
ANOVA program RoT 
Calculate ANOVA 
table entries A ss. (1) 
Note: the numbers in 
parentheses in the R/S r-1 (2) 
display column 
correspond to the R/S MSS 63) 
numbered blocks in 
the ANOVA table R/S 
(table 4-4) 
R/S 
R/S 
R/S 
R/S 
R/S 
R/S Prob-value 
(10) 
TABLE B-13.1 
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ANOVA USER INSTRUCTIONS 


Confidence Intervals fo} 
Contrasts 


10 Initialize 


ll Enter contrast data 


Repeat for each row 1 


eZ Enter F percentile with 
degrees of freedom r-l, 
ein) 
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CATEGORY I - ADJUSTED STRENGTH SCORES 
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CATEGORY I - ADJUSTED ENDURANCE SCORES 
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CATEGORY II - ADJUSTED STRENGTH SCORES 
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CATEGORY II - ADJUSTED ENDURANCE SCORES 
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CATEGORY III - ADJUSTED STRENGTH SCORES 


BENCH PRESS + CURL 
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CATEGORY III - ADJUSTED ENDURANCE SCORES 
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CATECORY IV - ADJUSTED STRENGTH SCORES 


BENCH PRESS + CURL 
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CATEGORY IV - ADJUSTED ENDURANCE SCORES 
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CATECORY V - ADJUSTED STRENGTH SCORES 
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CATEGORY V - ADJUSTED ENDURANCE SCORES 
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CATEGORY VI - ADJUSTED STRENGTH SCORES 
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CATEGORY VI - ADJUSTED ENDURANCE SCORES 
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